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Abstract. We construct eigenspectra from the DRl quasars in the 
SDSS using the Karhunen-Loeve (KL) transform (or Principal Compo- 
nent Analysis, PCA) in different redshift and luminosity bins. We find 
that the quasar spectra can be classified, by the first two eigenspectra, into 
a continuous sequence in the variation of the spectral slope. We also find a 
dependence on redshift and luminosity in the eigencoefficients. The domi- 
nant redshift effect is the evolution of the blended Fe II emission (optical) 
and the Balmer continuum (the "small-bump", Xrest ~ 2000 - 4000A), 
while the luminosity effect is related to the Baldwin effect. Correlations 
among several major broad emission lines are found, including the well- 
known "Eigenvector-1" . 



1. Why Study Quasar Eigenspectra? 

Spectroscopic observations from the Sloan Digital Sky Survey (SDSS) have the 
advantage of large numbers of quasars (QSOs) and a large redshift range that 
provides an unique opportunity to study their intrinsic properties, sample vari- 
ation, and spectral classification in great detail. A powerful method to address 
these questions is the KL transform. The idea is to derive from the observed 
spectra a lower dimensional set of eigenspectra (Connolly et al. 1995), in which 
the essential spectral properties are represented. In the following we outline 
some results from applying this method to the SDSS quasars. 



2. Data and KL Transform 

The sample we use is the Sloan Digital Sky Survey (SDSS; York et al. 2000) 
First Data Release (DRl; Abazajian et al. 2003) quasar catalog (Schneider et 
al. 2003), consisting of 16,713 quasars. The redshifts range from 0.08 to 5.41 and 
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Figure 1. Observed quasar spectra for different values of the first and 
second eigencoefficients (in a particular redshift and luminosity bin). 



the K-corrected i-band absolute magnitudes (Mi) from -30 to -22. We divide 
the whole sample into several redshift and luminosity bins with Az ranging from 
~ 0.45 to 1.8 and AMj = 2. We consider wavelength ranges for all bins within a 
restframe 900 A to 8000 A. The spectra are shifted to their restframes, and bad 
data are approximated using the gap-repairing formalism (Connolly & Szalay 
1999). The details of the KL transform are described in, for example, Connolly 
et al. (1995). Previous work applying the KL transform on quasar spectra or 
observed parameters include Francis et al. (1992), Boroson & Green (1992), and 
Shang et al. (2003). In the following, the flux densities are expressed in restframe 
values and the wavelengths are in vacuum. 



Eigenspectra of SDSS DRl QSOs 



3 



0.1 



0.08 



0.06 



0.04 



0.02 



-0.02 



-0.04 



4000 



1st eigenspectrum (= mean spectrum) 
2nd eigenspectrum 




4500 5000 
restframe wavelength (A) 



5500 



Figure 2. 
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Eigenspectra of the SDSS DRl quasars around the H/3 



Results 



We find that the quasars can be classified into a continuous variation of the 
spectral slope. Figure 1 shows the observed quasar spectra as a function of the 
first two eigencoefficients, ai and 02 (in our notation the first mode is the mean 
spectrum), in the redshift range 1.16 — 2.06 and the luminosity range —28 to 
—26. The actual values of the eigencoefficients are shown on top of each sub- 
figure. Along the sequence with decreasing 02 values, the quasar continua are 
progressively bluer. This implies that the linear-combination of the first 2 modes 
changes the spectral slope. In this regard, this is similar to the galaxy spectral 
classification by the KL approach (Connolly et al. 1995). 

Figure 2 shows the first two eigenspectra in the local region around H fj 
(A^est = 4000 - 55OOA), constructed with more than 2,500 SDSS QSOs in the 
redshift range 0.08 — 0.67. An anti-correlation is seen between Fe II (optical) 
and [O III]A5008, which is also known as the "Eigenvector-1" (Boroson &; Green, 
1992). It is generally believed that the [O III]A5008 luminosity is an isotropic 
property, as such this anti-correlation is not driven by the observed orientation 
of the quasar. 

In the local region around Ly a {Xrest = 1150 — 2000A), the eigenspectra 
from about 1300 wavelength-selected SDSS QSOs show agreement with those 
from Francis et al. (1992). The second eigenspectrum shows low- velocity (so- 
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called line-core) components of broad emission lines, and the third mode gives 
the continuum slope variation. 

We find a dependence of the coefficient from the second eigenspectrum 
on both redshift and luminosity (not shown). The dominant redshift effect is 
found to be the evolution of the blended Fc II emission lines (optical) and the 
Balmer continuum (the "small-bump" , \rest ~ 2000 — 4000A) , as such it is more 
prominent in the lower-redshift quasars and less obvious or absent in the higher- 
redshift ones. This result suggests that the Iron abundances build up as quasar 
evolves, in agreement with quasars in the Large, Bright Quasar Survey (Green 
et al. 2001) and in another independent study (Kuhn et al. 2001). On the other 
hand, the dominant luminosity effect in the quasar spectra is found to be the 
Baldwin effect — the decrease of line equivalent width with luminosity, which 
is detected in Ly a, Si IV+0 IV], C IV, He II, C III] and Mg II emission lines. 
The effect in N V seems to be redshift dependent, so that an opposite sign (an 
increase of line equivalent width with luminosity) is seen in the lower-redshift 
quasars. 

In summary, we show that the quasar eigenspcctra are invaluable to us in 
probing the intrinsic properties of quasars in a large data set. In addition to the 
redshift and luminosity dependent sets of eigenspcctra, we also construct a single 
set of cigcnspectra of the full DRl catalog ranging from restframe wavelengths 
900A to 8OOOA in order to study the correlations between the UV and the optical 
spectral regions. Interested readers can find detail discussions in properties of 
these cigcnspectra and future directions in Yip et al. (2003). 
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